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AND METHOD FOR ESTIMATING CHARGE RATE OF 
SECONDARY CELL 

BACKGROUND OF THE INVENTION: 

Field of the invention 

[0001] The present invention relates to apparatus 
and method for estimating a charge rate (abbreviated 
10 as SOC) of a secondary cell. 

Description of the related art 

[0002] Japanese Patent Application First 

15 Publications No. 2000-323183 published on November 24, 
2000 No. 2000-268886 published on September 29, 2000, 
and a Japanese Paper titled Estimation of Open 
Voltage and Residual Values for Pb Battery by 
Adaptive Digital Filter announced by a Japanese 

20 Electrical Engineering Society (T.IEEE Japan), Volume 
112-C, No. 4, published on 1992 exemplify previously 
proposed SOC estimating apparatus for the secondary 
cell. That is to say, since the charge rate (or 
called State Of Charge, i.e., SOC) of the secondary 

25 cell has a correlation to an open-circuit voltage Vq 
(cell terminal voltage when its power supply of the 
cell is turned off, also called electromotive force 
or open voltage), the charge rate can be estimated 
when open voltage Vq is obtained. However, a 

30 considerable time is needed until the terminal 
voltage is stabilized after the power supply is 
turned off ( charge-and-discharge is ended) . Hence, a 
predetermined time duration is needed from a time at 
which the charge-and-discharge is ended to determine 
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an accurate open-circuit voltage Vq . Therefore, since 
immediately after or during the charge/discharge time 
or charge-and-di scharge , it is impossible to 
determine an accurate open-circuit voltage and the 
5 charge rate cannot be obtained using the above- 
described method. Nevertheless, to determine the 
open-circuit voltage Vq, the open-circuit voltage Vq 
is estimated using a method disclosed in the above- 
described Japanese Patent Application First 
10 Publication No. 2000-323183. 

SUMMARY OF THE INVENTION: 
[0003] However, in the above-described method 

disclosed in the Japanese Patent Application 

15 Publication No. 2000-323183, open-circuit voltage Vq 
is calculated from a non-recursive (non-regression 
type) cell model (a model whose output value is 
determined only from a present value and past value 
of an input value) whose characteristic is wholly 

20 different from a physical characteristic of the cell 
for. which an adaptive digital filter (sequential type 
model parameter identification algorithm) is used. 
The charge rate (SOC) is used from this value. Hence, 
in a case where this method is applied to the actual 

25 cell characteristic (input: current, output: voltage), 
according to the cell characteristic, an estimation 
calculation is wholly converged or does not converge 
to a real value. Hence, it is difficult to estimate 
the charge rate (SOC) accurately. 

30 [0004] It is, hence, an object of the present 

invention to provide apparatus and method for 
estimating accurately the charge rate (SOC) for the 
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secondary cell and accurately estimating other 
parameters related to the charge rate (SOC) . 
[0005] According to one aspect of the present 

invention, there is provided a charge rate estimating 
5 apparatus for a secondary cell, comprising: a current 
detecting section capable of measuring a current 
flowing through the secondary cell; a terminal 
voltage detecting section capable of measuring a 
voltage across terminals of the secondary cell; a 

10 parameter estimating section that calculates an 

adaptive digital filtering using a cell model in a 
continuous time series shown in an equation (1) and 
estimates all parameters at one time, the parameters 
corresponding to an open-circuit voltage Vo, which is 

15 an offset term of the equation, (1) and coefficients 
of A(s), B(s), and C(s), which are transient terms; 
and a charge rate estimating section that estimates 
the charge rate from a previously derived 
relationship between an open-circuit voltage and a 

20 charge rate of the secondary cell and the open- 
circuit voltage Vo , 

V = ^^^•/+^— •Fn (l)f wherein s denotes a 

A(s) C(5) ' 

Laplace transform operator, A(s), B(s), and C(s) 

25 denote poly-nominal functions of s. 

[0006] According to another aspect of the present 

invention, there is provided a charge rate 
estimating method for a secondary cell, comprising: 
measuring a current flowing through the secondary 

30 cell; measuring a voltage across terminals of the 
secondary cell; calculating an adaptive digital 
filtering using a cell model in a continuous time 
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series shown in an equation (1); estimating all 
parameters at one time, the parameters corresponding 
to an open-circuit voltage Vq, which is an offset 
term of the equation (1), and coefficients of A(s), 
5 B(s), and C(s), which are transient terms; and 

estimating the charge rate from a previously derived 
relationship between an open-circuit voltage and a 
charge rate of secondary cell, and the open-circuit 
voltage Vo, 

10 V = ^^^^ • /+ — ^ — •Fn (l)f wherein s denotes a 

A(s) C(s) ' 

Laplace transform operator, A(s), B(s), and C(s) 
denote poly-nominal functions of s. 

[0007] According to a still another object of the 

present invention, there is provided a charge rate 

15 estimating method for a secondary cell, comprising: 
measuring a current I(k) flowing through the 
secondary cell; measuring a terminal voltage V(k) 
across the secondary cell; storing the terminal 
voltage V(k) when a current is zeroed as an initial 

20 value of the terminal voltage AV(k) = V(k) - V_ini; 

determining instantaneous current values Io(k), Ii(k), 
and l3(k) and instantaneous terminal voltages Vi(k), 
V2(k), and V3(k) from an equation (19), 

1 

lo = •// 

25 II =77^*/. Vx = -Zi^^^Vr 

12 = V2 = •V , 

13 = •/, V3 = •K, and 
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1 



1 



(19) , wherein pi denotes a 



constant determining a responsive characteristic of 
Gi(s); substituting the instantaneous current values 
Io(k), Ii(k), Iz (k), and l3(k) and the instantaneous 
terminal voltages Vi(k), V2(k), and V3(k) into an 
equation ( 18 ) , 

^k) 



Y{k) 

e (k) 

P (k) 
10 Xi{k) = 



e(k -1) - y ( k) •? (k-l) 'CO (k) • [co^ (k) (k -1> -y(k)] 



I A,ik)*Pik-l)»(oik)*a>''(k)*P(k-iy 

A,(A:)1 l + Z3ik)»co''(k)»P(k-V)»a)ik) 



^ trace {p\k) } ^ ^ trace 



{P\k) } 



Yu 



trace 



{P\k)} 



Yu 



trace 



{P\k) } 



Yl 



trace{p\k) } . trace{p' jk) } ^ 



Yl 



Yl 



15 -(18), 

wherein 9(k) denotes a parameter estimated value at a 
time point of k (k = 0, 1, 2, 3 ---), A-i, A,3(k), yu , 
and yL denote initial set value, b < A,i < 1 , 0 < A,3(k) 
< 00. P(0) is a sufficiently large value, G(0) 

20 provides an initial value which is non-zero but very 
sufficiently small value, trace{P} means a trace of 
matrix P, wherein y{k) = Vi(k) 

co^(k) = [V3(k) V2(k) l3(k) l2(k) Ii(k) 

Io(k)] 
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10 



15 



e (k) = 



(20); 



--a{k) 

c(k) 
d{k) 
e{k) 

J 

substituting a, b, c, d, e, and f in the parameter 
estimated value 9(k) into and equation (22) to 
calculate Vq which is an alternate of Vq which 
corresponds to a variation AVo(k) of the open-circuit 

voltage estimated value from a time at which the 
estimated calculation start is carried out; 

Vo' = (^1 ^^'^^\ vq = a • Ve + b • Vs +V4 -c • le - d • I5 - 

e • I4 (22); and calculating an open-circuit 

voltage estimated value Vo(k) according the variation 
AVo(k) of the open-circuit voltage estimated value 
and the terminal voltage initial value V_ini, 
[0008] This summary of the invention does not 

necessarily describe all necessary features so that 
the invention may also be a sub-combination of these 
described features. 



BRIEF DESCRIPTION OF THE DRAWINGS: 
[0009] Fig. 1 is a functional block diagram of an 

20 apparatus for estimating a charge rate (SOC) of a 

secondary cell in a preferred embodiment according to 
the present invention. 

[0010] Fig. 2 is a specific circuit block diagram 
of the apparatus for estimating the charge rate (SOC) 
25 of the secondary cell in the preferred embodiment 
according to the present invention. 
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[0011] Fig. 3 is a model view representing an 
equivalent circuit model of the secondary cell. 
[0012] Fig. 4 is a correlation map representing a 
correlation between an open-circuit voltage and a 
5 charge rate (SOC) . 

[0013] Fig. 5 is an operational flowchart for 
explaining an operation of a microcomputer of a 
battery controller of the charge rate estimating 
apparatus in the first preferred embodiment shown in 
10 Fig. 1. 

[0014] Figs. 6A, 6B, 6C, 6D, 6E, 6F, 5G, 6H and 61 

are characteristic graphs representing results of 
simulations of current^ voltages, and various 
parameters in a case of the charge rate estimating 
15 apparatus in the embodiment shown in Fig. 1 

DETAILED DESCRIPTION OF THE PREFERRED EMBODIMENTS: 
[0015] Reference will hereinafter be made to the 
drawings in order to facilitate a better 

20 understanding of the present invention. 

[0016] Fig. 1 shows a functional block diagram of 

charge rate estimating apparatus in a first preferred 
embodiment according to the present invention. In Fig. 
1, a reference numeral 1 denotes a parameter 

25 estimating section based on a cell model with an 

open-circuit voltage Vq (k) as an offset term. In 
addition, a reference numeral 2 denotes a open- 
circuit voltage calculating section to calculate 
open-circuit voltage Vq (k), and a reference numeral 

30 3 denotes a charge rate estimating section that 
calculate the charge rate from the open-circuit 
voltage. In addition, a reference numeral 4 denotes 
a current I measuring block to detect current I(k) 
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which is charged and discharged into and from the 
cell, and a reference numeral 5 denotes s terminal 
voltage of the cell to measure the terminal voltage 
V{k) . 

5 [0017] Fig. 2 shows a block diagram representing a 

specific structure of the charge rate estimating 
apparatus in the first embodiment. In this embodiment, 
a load such as a motor is driven with the secondary 
cell and the charge rate estimating apparatus is 

10 mounted in a system to charge the secondary cell with 
a regenerative power of the motor (load) . In Fig. 2, 
a reference numeral 10 denotes a secondary cell 
(simply called, a cell), a reference numeral 20 
denotes a load such as a DC motor, a reference 

15 numeral 30 denotes a battery controller (electronic 
control unit) to estimate the charge rate (charge 
state) of the cell having a microcomputer including a 
ROM (Read Only Memory) , a RAM (Random Access Memory) , 
a CPU (central Processing Unit), and Input/Output 

20 Interface and other electronic circuits. A reference 
numeral 40 denotes a current meter to detect a 
current which is charged into or discharged from the 
cell, a reference numeral 50 denotes a voltage meter 
to detect the terminal voltage of the cell, a 

25 reference numeral 60 denotes a temperature meter to 
detect a temperature of the cell. These meters are 
connected to battery controller 30. Battery 
controller 30 corresponds to parts of parameter 
estimating section 1, an open-circuit voltage Vq (k) 

30 and a charge rate estimating section 3. Current 
meter 40 corresponds to current I(k) measuring 
section and voltage meter 50 correspond to terminal 
voltage V(k) measuring section 5. 
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[0018] First, a " cell model " used in the first 

embodiment will be described below. Fig. 3 is an 
equivalent circuit representing an equivalent circuit 
model of the secondary cell. The equivalent circuit 
5 model of the secondary cell can be represented by the 
following equation (7) (= equation (6)). 

V = ^ + Vo (7). 

In equation (7), a model input is a 

10 current I [A] (a positive value represents a charge 

and a negative value represents a discharge), a model 
output is a terminal voltage V [V] , an open-circuit 
voltage is Vo, K denotes an internal resistance, Ti 
through T3 denote time constants (Ti T2 ^ T3, Ti« T3) 

15 and s denotes a Laplace transform operator. 

[0019] In this model based on equation of (7) is a 
reduction model (first order) in which a positive 
pole and a negative pole are not specially separated 
from each other. However, it is possible to represent 

20 a charge-discharge characteristic of an actual cell 

relatively easily. Equation (7), in equation (1) of V 

= B(s)/A(s)»I + l/C(s)*Vo (1), A{s) = Ti • s + 1, 

B(s) = K • (T2 • s + 1), C(s) - T3 • s + 1. 

[0020] Hereinafter, a deviation from the cell model 

25 based on equation (7) to an adaptive digital filter 
will first be described below. Open-circuit voltage 
Vo can be described by an equation (8), supposing 
that a value of a current I multiplied with a 
variable efficiency of A is integrated from a certain 

30 initial state. 

That is to say, Vo = — •/ (8) . 

s 
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It is noted that equation (8) corresponds to 
a replacement of h recited in equation (2), viz., Vq 
= h/s • I with efficiency of A. 

[0021] If equation (8) is substituted into equation 
5 (7), equation (9) is resulted. 



Vo = — + • — •/ (9). 

r, •5 + 1 T^^s + l s 

Equation (9) corresponds to equation (3) (V 



\^A(s) C(s) s) s^A{s)•C{s) 
10 A(s), B(s), and C(s) in equation (3), the following 
equations are substituted into equation (9) in the 
same way as the case of equation (7) . 
A ( s ) = Ti • s + 1 , 
B (s) = K • (T2*s + 1) 
15 C(s) = T3*s -I- 1. In other words, equation 

(3) is a generalized equation and this application to 
a first order model is equation (9) . If equation (9) 
is arranged, an equation of (10) is given. 
S •(Ti»s + l)(T3»s + l)*V = KMT2*s + 1) (Ta^s + l)s»I + 
20 AMTi^s + 1)»I 

{Ti* Ta* s^ + (Ti + T3)»s^ + s} • V = {K#T2»T3«s^ + K • 
{T2 + T3)*s^ + (K + A • Ti)*s + A} • I 

{a»s^ + h^s^ + s)*V = {c^s^ + d*s^ -h e • s + f)*I 

(10). It is noted that, in the last equation of 
25 equation (10), parameters are rewritten as follows: 
a = Ti#T3, b = Ti + T3, c = K • T2 • T3, d = K • (T2 + 
T3) , e = K + A • Ti, and f = A (11)- 
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If a stable low pass filter Gi(s) is introduced into 
both sides of equation (10) and arranged, the 
following equation (12) is given. 

— ^- — (a^s^ + b*s^ + s)*V = — ^ — (c»s^ + d«s^ + e^s + f)*I 

5 (12). In details, in equation (10), on the 

contrary of equation (7), if Ti • s + 1 = A(s), K* 
(T2»s + 1) = B(s), and Ta-s + 1 = C(s) are substituted 
into equation (10), this is given as: s*A ( s ) ( s ) = 
B ( s ) •C ( s ) •s* I + A»A(s)»I. This is rearranged as 
10 follows: s#A ( s ) #0 ( s ) •V = [B(s)«C(s)»s [['I]] + 

A«A(s)]-I (12)', If, the low pass filter (LPF), 

Gi(s) is introduced into both sides of equation (12) 
an equation (4) is given. 

15 That is to say, — ^•V = ',\',^ ^ ' •! 

G,is) s^A{s)•Cis) 

(4). It is noted that s denotes the Laplace 
transform operator, A(s), B(s), and C(s) denote a 
poly-nominal function of s, h denotes a variable, and 
1/Gi(s) denotes a transfer function having a low pass 

20 filter characteristic. That is to say, equation (4) 
is the generalized function, equation (12) is the 
application of equation (4) to the first order model. 
[0022] Current I and terminal voltage V which can 

actually be measured are processed by means of a low 

25 pass filter (LPF) and a band pass filter (EPF) are 
defined in the following equations (13), provided 
that pi denotes a constant to determine a responsive 
characteristic of Gl(s) and is determined according 
to a designer's desire. 
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lo = 



Ii = 



I2 = 



I3 = 



1 



1 

G,(5) 

s 



•I, 



•I, 



1 



G,(5) (F,*s + iy 



Vi = 



V, = 



G,(5) 



V3 = 



G,(s) 



•V , 



(13) 



10 



15 



If equation (12) is rewritten using the variables 
shown in equations (13), equations (14) are 
represented and, if deformed, the following equation 
(15) is given . 

a»V3 + b»V2 + Vi = c»l3 + d»l2 + e«Ii + 

f • lo 

Vi = -a*V3 - b»V2 + c»l3 + d*l2 + e • Ii + 
f»Io (14 ) . 



-a 
-b 

c 
d 
e 

f . 



Vi = [V3 V2 I3 I2 Ii lo] • [ [ = ] ] 



(15) . 



20 



Equation (15) is a product-sum equation of measurable 
values and unknown parameters. Hence, a standard 
(general) type (equation (16)) of the adaptive 
digital filter is coincident with equation (15). 
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It is noted that (o*^ means a transposed vector in 
which a row and column of a vector (o are mutually 
exchanged . 

y = co^*0 (16) . It is noted that y, 



(0 



and 0 can be expressed in the following equation (17) 
in equation (16) described above. 



Y = Vi, co^ = [V 



V2 



l2 



lo] . e = 



(17) 



- a 
-b 
c 
d 
e 

/ . 

Hence, if a signal filter processed for current I and 
terminal voltage V is used in a digital filter 
10 process calculation, unknown parameter vector 9 can 
be estimated. 

[0023] In this embodiment, " a both-limitat ion 

trace gain method is used which improves a logical 
demerit of a simple " an adaptive digital filter by 

15 means of a least square method " such that once the 
estimated value is converged, an accurate estimation 
cannot be made any more even if the parameters are 
changed. A parameter estimating algorithm to estimate 
unknown parameter vector 0 with equation (16) as a 

20 prerequisite is as shown in an equation (18). It is 
noted that the parameter estimated value at a time 
point of k is 0(k) . 



y (k) = 



\^X^{k)^co^ (k) • P(k - 1) • (Oik) 



25 0(k) = 0(k -1) - Y ( k) ( k-1 ) •© ( k) • [co^ ( k) •© ( k -1> ~y(k)] 
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Xi(k) = 



r 



v. 



trace {p'(k) } trace{p'ik) } 

ru Yu 

trace{p\k) } ^^^^ trace{p\k) } 
trace{p\k) } . trace{p\k) } ^ 



(18) . 

In equations (18), Xi, X3(k), Yu, and Yl 
denote initial, set value, b < A,i < 1 , 0 < X3(k) < oo . 
P(0) is a sufficiently large value, 9(0) provides an 

10 initial value which is non-zero but very sufficiently 
small value. In addition, trace{P} means a trace of 
matrix P. As described above, the derivation of the 
adaptive digital filter from cell model. 
[0024] Fig. 5 shows an operational flowchart 

15 carrying out the microcomputer of battery controller 
30. A routine shown. in 5 is carried out for each 
constant period of time Tq. For example, I(k) is the 
present value and I(k -1) means a one previous value 
of I(k). At a step SIO, battery controller 30 

20 measures current I(k) and I(k -1) means one previous 
value of I(k). At a step S20, battery controller 30 
carries out a turn on-and-off determination of an 
interrupt relay of the secondary cell. That is to 
say, battery controller 30 performs the on-and-off 

25 control of the interr.upt relay of the secondary cell- 
When a relay is turned off (current I = 0), the 
routine goes to a step S30. During the engagement of 
the relay, the routine goes to a step S40. At step 
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S30, when the relay is engaged, the routine goes to a 
step S540. At step S530, battery controller 30 
serves to store terminal voltage V(k) to as an 
initial value of the terminal voltage V_ini . At a 
5 step S40, battery controller 30 calculates a 

differential value AV(k) of the terminal voltage. 
AV(k) = V(k) - V_ini . This is because the initial 
value of the estimation parameter in the adaptive 
digital filter is 0 so that the estimation parameter 

10 does not converge during the estimation calculation 
start time. Thus, all of inputs are zeroed. During 
the input being all zeroed. During the relay 
interruption, step S30 have been passed and the 
estimation parameters are remains initial state since 

15 1=0 and the estimation parameter remains alive. 

[0025] At step S50, a low pass filtering or band 

pass filtering are carried out the current I(k) and 
terminal voltage difference value AV ( k ) on the basis 
of equation (13). Io(k) through l3(k) and Vi(k) 

20 through V3(k) are calculated from equation (19) . In 
this case, in order to improve an estimation accuracy 
of the parameter estimation algorithm of equation 
(18), a responsive characteristic of low pass filter 
Gi(s) is set to be slow so as to reduce observation 

25 noises. However, if the characteristic is quicker 

than a response characteristic of the secondary cell 
(a rough value of time constant Ti is known), each 
parameter of the electric cell model cannot 
accurately be estimated. It is noted that pi recited 

30 in equation (19) denotes a constant determined 

according to the responsive characteristic of Gi(s). 
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[0026] 



G,(5) 



G,(5) 



10 = 

11 = 

12 = 

13 = 
1 _ 1 



•I, 



•I, 



Vi = 



V2 = 



V3 = 



„3 



•V , 



G,(5) 



• V , and 



(19) 



At a step S60, Io(k) through l3(k) calculated at step 
S50 and Vi(lc) through V3(k) are substituted into 
equation (18) - Then^ the parameter estimation 
10 algorithm in the adaptive digital filter^ viz., 
equation (18) is executed to calculate parameter 
estimated value 0(k). y(k), co^(k), and 9(k) are shown 
in equation (20) . 

y (k) = Vi (k) 

15 co^(k) = [V3(k) V2(k) l3(k) l2(k) Ii(k) 

Io(k) ] 



e(k) = 



c{k) 
d{k) 
e{k) 



(20) . 



At a step S70, a through e of parameter estimated 
20 value 9(k) calculated at step S50 are substituted 

into the following equation (22) in which the above- 
described cell model equation (7) is deformed to 
calculate Vq ' which is an alternative to open-circuit 
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voltage Vo. Since the variation in open-circuit 
voltage Vo is smooth, Vo ' can be used alternatively. 
It is noted that the derivation herein is a variation 
AVo(k) of the open-circuit voltage from the estimated 
5 calculation start time. 

[0027] It is noted that an equation of [l/Cl(s)]I 

in equation (21) is replaced with an equation (24) 
corresponds to equation (22) . It is also noted that, 
in the derivation of equation (22), K in equation 
10 (21) is strictly different from e in equation (21) . 
However, since, physically, K » A • Ti, e is 
approximated to K (e = K) . Then, each coefficient a 
through e in equation (22) is the contents shown in 
equation (23). 

15 •Fq = V - ^ — r^^i 

(Ti»s + !)• Vq =(Ti»s + l)(T3#s + 1)V - 

K^ (T2»S + 1 ) (T3*S + 1 ) •I 

(Ti*s + 1)* Vo = {Ti#T3«s^ + (Ti + T3)»s + 

1}#V 

20 - {K»T2»T3«S^ + K«(T2 +T3)»s + K}»I 

(21) . 

V'o = (^1 • + 1) ^ ^ a^Ve + b»V5 +V4 -c^Ie - d#l5 - e» I4 - 

G2{s) 

(22) . 



25 

[0028] It is noted that a = Ti • T3. b = Ti + T3, c 

= K •(T2 + T3}, d = K •(T2 + T3), e = K + A#Ti = K-- 

- (23) . 
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[0029] l4 = 77^^^. V4 = T^^^' 

I5 = -^•I , V5 = -^•V , 
1 1 1 



Is = ^•f^ and Ve = -^-V (24). 

5 [0030] P2 recited in equations (24) denote a 

constant to determine a responsive characteristic of 
G2(s) . Ti of the cell parameter is known to be 
several seconds. Hence, T'l in equation (24) is set 
to be approximated value to Ti- Thereby, since (Ti* s 

10 + 1 ) which remains in a numerator of equation (22) 
can be compensated, the estimation accuracy of open- 
circuit voltage Vq can be improved. It is noted that 
equation (21) corresponds to equation (5) . That is to 
say, equation (21) is derived from (Ti • s +!)• Vo = 

15 (Ti • S + 1) (T3 • s + 1)^V - K*(T2* s + 1) (Ts* S +1)* 

(T3 • s + 1)* I. If the following three equations are 
substituted into the above-described deformation of 
equation (21). Ti • s + 1 = A(s), K • (T2 • s + 1) 
B(s), and T3 • s + 1 = C(s). That is to say, A{s)^ Vq 

20 = A(s)* C(s) • V - B(s)^ C(s)^ I. If this is 

rearranged, this results in Vq = C(s) • V - B(s) • 
C(s) • I/A(s), Vo = C(s)#[V - B(s) • I/A(s)] If low 
pass filter G2(s) is introduced into both sides of 
this equation, this results in equation (5) . In 

25 details, equation (5) is a generalization equation 
and the application of equation (5) to the first 
order model is equation (2) . 
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[0031] At a Step S80, battery controller 30 adds 

the open-circuit voltage initial value, i.e., 
terminal voltage initial value V_ini to a variation 
AVo(k) of open-circuit voltage Vq so as to obtain 
5 open-circuit voltage estimated value Vo(k) from the 
following equation (25) . 

Vo (k) = A Vo (k) + V_ini (25). 

[0032] At a step S90, battery controller 30 

calculates the charge rate SOC(k) from open-circuit 

10 voltage Vq (k) calculated at step S80 using a 

correlation map of the open-circuit voltage versus 
the charge rate as shown in Fig. 4. It is noted that, 
in Fig. 4, Vl denotes the open-circuit voltage 
corresponding to SOC = 0 % and Vh denotes the open- 

15 circuit voltage corresponding to SOC = 100 %. • At a 

step SlOO, battery controller 30 stores the necessary 
numerical values needed in the subsequent calculation 
and the present routine is ended. As described above, 
an operation of the apparatus for estimating the 

20 charge rate of the secondary cell has been described. 
[0033] (1) As described above, a relationship 

from among current I of the secondary cell and 
terminal voltage V thereof, and the open-circuit 
voltage Vq is structured in transfer function that as 

25 in the general equation (1), that in the preferred 
embodiment, equation (7) (= equation (6). Hence, it 
is made possible to apply an adaptive digital filter 
such as a least square method (well known estimation 
algorithm) . Consequently, it becomes possible to 

30 estimate parameters in equations (viz., open-circuit 
voltage Vq which is an offset term and poly-nominal 
equations A(s), B(s), and C(s)) in a form of a batch 
processing. These parameters are largely affected by 
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the charge rate, a surrounding temperature, and a 
deterioration and varied instantaneously. It is 
possible to sequentially estimate the adaptive 
digital filter with good accuracy. Then, if a unique 
5 correlation between the open-circuit voltage Vo and 
the charge rate as shown in Fig. 4 are stored, the 
estimated open-circuit voltage can be converted to 
the charge rate. Hence, it is possible to 
sequentially estimate the charge rate in the same way 

10 as the parameters described above. 

[0034] (2) In a case where the equation (1) which 

is the relationship equation of current I and 
terminal voltage V of the secondary cell is 
approximated to equation (4), the equation such that 

15 no offset term is included (viz., the open-circuit 

voltage Vo ) / a product-and-addi t ion equation between 
a measurable current I which is filter processed and 
a terminal voltage V which is filter processed and 
unknown parameter (coefficient parameters of poly- 

20 nominal equations A(s), B(s), and C(s) and h) is 
obtained. A normally available adaptive digital 
filter (the least mean, square method and well known 
parameter estimation algorithm) can directly be 
applied in a continuous time series. 

25 [0035] As a result of this, the unknown parameters 

can be estimated in the batch processing manner and 
the estimated parameter h is substituted into 
equation (2), the estimated value of open-circuit 
voltage Vo can easily be calculated. All of these 

30 parameters are varied instantaneously, the adaptive 
digital filter can serve to estimate the charge rate 
at any time with a high accuracy. Since a constant 
relationship between open-circuit voltage Vo and the 
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charge rate (SOC) is established as shown in Fig. 4, 
if this relationship is previously stored, the charge 
rate (SOC) can be estimated from the estimated value 
of open-circuit voltage Vq - 
5 [0036] Figs- 6A through 61 integrally shows 

signal timing charts with current I and terminal 
voltage V inputted into adaptive digital filter and 
representing results of simulation graphs when each 
parameter is estimated. As far as a time constant of 

10 a first order delay in equation (6) is concerned, Ti 
< To. Since all parameters a through f (refer to 
equation (11)) are favorably estimated, the estimated 
value of open-circuit voltage Vq can be said to be 
well coincident with a real value. 

15 [0037] It is noted that, in Fig. 6C which 

indicates the open-circuit voltage, a reason that a 
right side second term of equation (6) is described 
is to indicate that the open-circuit voltage 
estimated value is coincident with a real value 

20 almost without delay in spite of the fact that a 
late term of time constant T3 is measured on the 
terminal voltage inputted into the adaptive filter. 
In details, since the parameter estimation with the 
cell model formatted adaptive digital filter in 

25 equation (5), all of parameters a through f can 

favorably be estimated and the estimated value of 
open-circuit voltage Vq is well coincident with a 
real value . 

[0038] (3) In addition, as described in item (2), 

30 in the structure in which the open-circuit voltage Vq 
is calculated from equation (2), the integration 
occurs before a value at which estimated value h is 
converged to the real value, its error cannot be 
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eliminated. However^ in the structure in which 
equation (5) in which the integration is not included, 
the error before the parameter estimated value is 
converged into the real value does not give an 
5 influence after the convergence. 

[0039] It will be appreciated that, in part of ® 

in Fig. 51, before estimated value f is converged 
into a real value, an erroneous estimation is carried 
out only at momentarily. In equation (2), this value 

10 is also integrated so that the error is not 

eliminated. However, the error is not eliminated 
since even this value is integrated. However, in the 
structure using equation (5), open-circuit voltage Vq 
is calculated from the equation in which the 

15 integration is not included. Hence, after the 

parameter estimated value is converged into the real 
time, this erroneous estimation portion is eliminated. 
[0040] (4) Furthermore, in a case where equation 
(6) is used in place of equation (1), a calculation 

20 time and program capacity can be suppressed to a 

minimum while having the above-described advantages. 
[0041] The entire contents of a Japanese Patent 

Application No. 2002-340803 (filed in Japan on 
November 25, 2002) are herein incorporated by 

25 reference. The scope of the invention is defined with 
reference to the following claims. 
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